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Abstract 

A total of 65 Bacillus thuringiensis (Bt) isolates were subjected to analysis of genetic relationship using fAFLP 
(fluorescent Fragment Length Polymorphism), in order to determine the genetic diversity within a group of Bt 
strains. 26 strains from different subspecies were identified as it follows: 9 kindly provided by the USDA (United 
States Department of Agriculture), 9 kindly provided by the Institute Pasteur and eight from Embrapa Maize and 
Sorghum Bt Collection, and 39 strains with no subspecies information also from Embrapa's Bt Collection. DNA 
sample was double digested with restriction enzymes fcoRI and Mse\, and the fragments were linked to adapters. 
Selective amplification reactions were performed using five primer combinations and the amplified fragments were 
separated by gel electrophoresis on an ABI377 sequencer. Genetic distances were obtained by the complement of 
the Jaccard coefficient and the groups were performed by the UPGMA method. Five primer combinations 
generated 495 scorable fragments and 483 were found to be polymorphic. Out of 26 subspecies, strains 344 and 
T09 (6. thuringiensis subsp. tolworthi) showed the highest similarity (15%), while isolates HD3 B. thuringiensis subsp 
finitlmus and T24 B. thuringiensis subsp neoleonensis were the most genetically distant (92%). B. thuringiensis isolates 
with no subspecies identification, found in samples from Goias State showed higher similarity forming a group with 
an average distance of 6%, and the closest subspecies to this group was B. thuringiensis subsp thuringiensis (HD2) 
with 52% of similarity. This similarity may be due to the fact that these organism exchange genetic material by 
conjugation, and it is relatively common to have evolutionary characteristics of their ancestors. 
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Background 

Bacillus thuringiensis (Bt) is a rod-shaped, Gram-positive 
soil bacterium characterized by its ability to produce crys- 
talline inclusions called "Cry proteins, Cry toxins or Bt 
toxins" during sporulation. The crystalline inclusions 
along with the spores have a high potential to control 
a great number of insect pests belonging to the order 
Lepidoptera, Diptera and Coleoptera (Vidyarthi et al. 
2002; Yamamoto and Dean 2000). B. thuringiensis can be 
found in different substrates such as soil, water, plant sur- 
faces, dead insects, grain dust, and stored grain (Glare and 
O'Callaghan 2000; Valicente and Barreto 2003), and it is a 
useful alternative or zment to synthetic chemical pesticide 
applications in commercial agriculture, horticulture, forest 
management, and mosquito control (De Maagd et al. 
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1999). B. thuringiensis also has high genetic variability and 
are widely distributed in nature (Vilas-Boas et al. 2002, 
Martin and Travers 1989, Hossain et al. 1997). The poly- 
merase chain reaction (PCR) has been used to identify 
specific cry genes and characterize B. thuringiensis strains 
with a high degree of success (Chal< et al. 1994; Ceron 
et al. 1994, 1995; Shangkuan et al. 2001; Lima et al. 2002; 
Valicente et al, 2010). 

The Amplified Fragment Length Polymorphism (AFLP) 
technique is based on the detection of genomic restriction 
fragments by PCR amplification, and can be used for 
DNAs of any origin or complexity (Vos et al. 1995). AFLP 
may be used as a tool for bacterial taxonomy and has 
shown utility in detecting molecular variability in very 
closely related bacterial strains (Jansen et al. 1997; Burke 
et al. 2004; Grady et al. 2001). The fAFLP (flurescent Amp- 
lified Fragment Length Polymorphism) use markers with 
fluorescent substances, the use of primers labeled with a 
fluorochrome, associated with automatic sequencers and a 
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large capacity for computational analysis (Ryu et al. 2005; 
Burke et al. 2004). 

Spodoptem frugiperda (Smith) is responsible for sig- 
nificant losses in maize production in Brazil. Its control 
is mainly achieved using chemical insecticides. However, 
biological control using B. thuringiensis and the produc- 
tion of genetically modified plants with insect resistance 
genes are excellent alternatives to control this insect 
(Valicente and Barreto 2003). 

The purpose of the present study was to use fAFLP for 
molecular analysis and characterization of B. thuringiensis 



*Mortality results from Valicente and Barreto 2003 and Valicente and Fonseca 2004. 



strains that were previously tested against S. frugiperdz, 
and to determine the degree of genetic variability and 
correlation between different genotypes of public strains 
of B. thuringiensis and strains sampled in different 
Brazilian regions. 

Materials and methods 

B. thuringiensis strains 

Of a total of 65 B. thuringiensis strains used in this 
research (Table 1), 26 strains from different subspecies 
were identified as it follows: 9 kindly provided by the 
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Table 1 Bacillus thuringiensis serovars used in this study, with mortality against fall armyworm, Spodoptera frugiperda 



N° 


Bt subspecies 


Strain 
identification 


Mortality (%)* 


Origin 


N° 


Strain 
identification 


Mortality (%)* 


Origin 


1 


alesti 


HD-4 


6.8 


USA 


34 


1109 N 


100 


Goias State/Brazil 


2 


alesti 


348B 


100 


Parana State/Brazil 


35 


1145B 


100 


Goias State/Brazil 


3 


aizawai 


T-07 


80.8 


France 


36 


1145C 


100 


Goias State/Brazil 


4 


aizawai 


HD-11 


7.8 


USA 


37 


1148 F 


100 


Goias State/Brazil 


5 


galleriae 


462A 


100 


Parana State/Brazil 


38 


1132E 


100 


Goias State/Brazil 


6 


galleriae 


474 


100 


Parana State/Brazil 


39 


1135B 


100 


Goias State/Brazil 


7 


galleriae 


348 L 


100 


Parana State/Brazil 


40 


1136B 


100 


Goias State/Brazil 


8 


galieriare 


HD-29 


12.8 


USA 


41 


1139 K 


100 


Goias State/Brazil 


9 


tolworthi 


344 


100 


Parana State/Brazil 


42 


BTLM 


100 


Goias State/Brazil 


10 


tolwonhi 


T-09 


100 


Frange 


43 


1644 


100 


Parana State/Brazil 


11 


tolworthi 


426 


100 


Ceara State/Brazil 


44 


1641 


100 


Parana State/Brazil 


12 


toiwortini 


461 A 


100 


Parana State/Brazil 


45 


701 A 


100 


Sao Paulo State/Brazil 


13 


darmstadiensis 


460 


100 


Parana State/Brazil 


46 


1658 


100 


Sao Paulo State/Brazil 


14 


darmstadiensis 


T-10 


77.9 


France 


47 


1646 


100 


Sao Paulo State/Brazil 


15 


entomocidus 


T-06 


9.8 


France 


48 


701 B 


100 


Sao Paulo State/Brazil 


16 


neoieonensis 


T-24 


1 7.9 


France 


49 


1648 


100 


Sao Paulo State/Brazil 


17 


mexicanensis 


T-27 


17 


France 


50 


1438 


100 


Sergipe State/Brazil 


18 


Japonensis 


T-23 


33.5 


France 


51 


1438H 


100 


Sergipe State/Brazil 


19 


Indiana 


T-16 


12.2 


France 


52 


1089E 


95.8 


Minas Gerais/State/Brazi 


20 


kurstaki 


HD-73 


2.7 


USA 


53 


1091H 


95.8 


Minas Gerais State/Brazi 


21 


kurstaki 


HD-1 


0 


USA 


54 


1604D 


95.2 


Amazonia State/Brazil 


22 


morrisoni 


HD-12 


28 


USA 


55 


1605 


95.2 


Amazonia State/Brazil 


23 


dendrolimus 


HD-7 


54 


USA 


56 


1657 


100 


Amazonia State/Brazil 


24 


finitlmus 


HD-3 


5.2 


USA 


57 


1656 


100 


Alagoas State/Brazil 


25 


tliuringlensis 


HD-2 


37.8 


USA 


58 


1503B 


100 


Santa Catarina State/Bra 


26 


Israelensis 


T-14 


0 


France 


59 


1626C 


97.6 


Maranhao State/Brazil 


27 




1119C 


100 


Goias State/Brazi 


60 


702 


97.6 


Mato Grosso State/Brazi 


28 




1124E 


100 


Goias State/Brazi 


61 


1354 


100 


Minas Gerais State/Brazi 


29 




1131'^ 


100 


Goias State/Brazi 


62 


1355SLO 


100 


Minas Gerais State/Brazi 


30 




1131C 


100 


Goias State/Brazi 


63 


1357E 


100 


Parana State/Brazil 


31 




1132^^ 


100 


Goias State/Brazi 


64 


566 


100 


Parana State/Brazil 


32 




1132C 


100 


Goias State/Brazi 


65 


1355LM 


100 


Minas Gerais State/Brazi 


33 




1138G 


100 


Goias State/Brazi 
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USDA (United States Department of Agriculture), 9 
kindly provided by the Institute Pasteur and eight from 
Embrapa Maize and Sorghum Bt Collection. Also, 39 
strains with no subspecies information also from Embra- 
pa's Bt Collection (Table 1). The words subspecies and ser- 
ovar were used interchangeably as given by Zeigler (1999). 
All strains have been previously tested against fall army- 
worm, S. frugiperda (Lepidoptera: Noctuidae J.E. Smith) 
(Valicente and Fonseca 2004; Valicente and Barreto 2003), 
and showed different toxicity activity. These strains are 
stored in glycerol in a freezer at - 20°C. 

Bacterial cultures 

B. thuringiensis isolates were inoculated into 5 ml LB 
medium enriched with salts and glucose (all expressed 
in g.r\ 0.002 g of FeS04, 0.02 g of ZnS04, 0.02 g of 
MnS04, 0.3 g of MgS04) and 2% glucose, with pH 
adjusted to 7.2. The shaking flasks were incubated at a 
stirrer speed of 250 rpm for 16 h at 28°C. Colonies were 
observed on a phase contrast microscope. 

DNA extraction 

Genomic DNA from B. thuringiensis strains was isolated 
and purified according to Shuhaimi et al. (2001) with 
modifications. A volume of 5 ml of an overnight culture 
was centrifuged at 4600 g for 5 min. The pellet was 
washed with sterile distilled water and resuspended in 
700 \A glucose/Tris/EDTA buffer (50 mM glucose, 25 mM 
Tris/HCl and 10 mM EDTA, pH 8). Lysozyme was added 
to the bacterial suspension at 20 mg/ml and incubated 
at 37°C for 30 min. The cells were additionally lysed by 
adding 15 [il 20% (w/v) SDS. Five |il of proteinase K at 
20 mg/ml was added and incubated at 65°C for 1 h to 
complete protein digestion. Five hundred \i\ of Phenol/ 
chloroform/isoamyl alcohol (25:24:1 by vol.) was added 
and centrifuged at 12.000 g for 1 min. Three hundred \i\ 
of the aqueous upper layer was then transferred into a 
sterile new tube and an equal volume of 3 M potassium 
acetate (pH 5.2), and a double volume of 2-propanol were 
added. The mixture was then centrifuged at 12.000 g for 
7 min. The supernatant was carefully discarded. The pellet 
was washed with 70% (v/v) ethanol by centrifuging at 
12.000 g for 5 min. The ethanol was discarded and the 
pellet containing the chromosomal DNA resuspended in 
150 |il of Tris-EDTA pH 8.0 (10 mM Tris-HCl, 1 mM 
EDTA) and stored at -20°C until further use. 

fAFLP analysis of DNA samples 

The fAFLP reactions were performed using to the Applied 
Biosystems AFLP™ microbial fingerprinting Kit with the 
genomic DNA (500 ng) being digested with 5U EcoRl and 
2,5 U oiMsel (Invitrogen", Carlsbad, CA, EUA) in 1,25 ^il 
of buffer lOx with final volume of 10 uL. After the diges- 
tions, were added 1 i^l of £coRI and Msel adapters 



(2 mM), 1.5 U of T4 DNA ligase buffer with T4 DNA 
ligase of (50 mM Tris-HCl, pH 7.8, 10 mM MgCl^ 
10 mM dithiothreitol, 1 mM ATP, 25 g/ml bovine serum 
albumin). The binding reactions were conducted at 20°C 
for 2 h, and amplicons produced were dissolved in 30 \A of 
TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8.0) and 
used as template in pre-selective amplification reactions. 
The binding sites of pre-selective primers are the se- 
quences of the adapters of restriction sites that provide a 
selection of restriction fragments obtained. The reactions 
were obtained using a 2 i^l of diluted DNA, 7.5 i^l of Core 
Mbc and 0.5 |il of EcoRl primers (1 |iM) and Msel (5 mM) 
were considered as pre-selective. Amplification conditions 
were: 72°C for 2 minutes, 20 cycles starting at 94°C for 
20 s, 56°C for 30 s and 72°C for 2 minutes. Products amp- 
lification were electrophoresed on 0.8% agarose gels and 
subsequently stained colored with ethidium bromide 
(0.5 mg/ml), and observations about the presence of bands 
at most of 500 bp. The pre-amplified products were 
diluted (1:15) with TE buffer and 1.5 \i\ were used in the 
pre-amplification reaction along with 7.5 i^l of buffer, 
nucleotides and AmpliTaq" DNA polymerase, a mixture 
called "Core Mbc, 0.5 \A of the Msel primer (5 \iM) and 
0.5 1^1 of EcoRl primer, labeled with fluorescence (1 |iM). 
Selective amplifications were performed in 20 i^l reactions 
with the following steps: 94°C for 30 s, 65°C for 30 s, and 
72°C for 1 min for 1 cycle and then lowering the annealing 
temperature by 1°C each cycle to 56°C, followed by an 
additional 20 cycles at a 56°C. Five primer combinations 
were used in the selective amplification {EcoRl -AC + 
Msel -CG; £coRI -G + Msel -CC; £coRI -C + Msel -G; 
£coRI -0 + Msel -CG; £coRI -C + Msel -CG). 

Sequencing eletrophoresis of polyacrylamide gel and 
data analysis 

Each sample was loaded onto a 5% Long-Ranger de- 
naturing gel which was run in an ABI PRISMTM 377 
DNA Sequencer (PE-Applied Biosystems) for 4 h at 
2,500 V. The electrophoresis images obtained from each 
reaction were analyzed using the GENESCAN program 
version 3.1 (Applied Biosystems). The Genotyper pro- 
gram version 2.5 (Applied Biosystems) was used to 

Table 2 AFLP primer combinations used for the 
fingerprinting of Bt serovars strains. 



Primer Number of Number of polymorphic 

combination fragments fragments 



* EcoRl- AC + Msel - CG 


77 


74 


* EcoRI-G + Msel- CC 


96 


89 


* EcoRl-C + Msel =G 


130 


129 


* EcoRl -0 + Msel -CG 


113 


112 


* EcoRl -C + Msel -CG 


79 


79 



*AFLP primer combinations. 
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check the presence (coded as 1) or absence (coded as 0) 
of polymorphic bands in the electropherograms, and 
produce a binary matrix and the GENES program ver- 
sion 2008 (Cruz 2001). The genetic relationships among 
strains were evaluated using a matrix of genetic dis- 
tances using the complement of the Jaccard similarity 
coefficient (CSJ), that does not consider negative simi- 
larities and the absence of the product (Cruz and Carneiro 
2006). From estimates of the dissimilarities, the lines were 



grouped using hierarchical UPGMA method (Unweighted 
Pair-Group Mean Average) with the test of bootstrap 
(1000 times) to evaluate the consistency of the group 
(Efron and Tibshirani 1993). 

Results and discussion 

In this study we have examined the genetic variation 
among isolates of B. thuringiensis. fAFLP markers 
generated excellent results in characterizing molecular 



_HD3 Btfinitimus - USA 
_J"24 Bt neoleonensis - Fran9a 
— T14 Bt israelensis — Franca 
-4ID4B(ates(i-USA 
_HD12 Bt morrisoni - USA 
-TIO Bt darmstadiensis - Franca 
-T27 Bt mexicanensis - Fran9a 
-T06 Bt entomocidus - Franca 
"T09 Bt tolworthi - Franca 

"344 Bt tolworthi - PR 
-^TLM - GO 
-1354 -MG 
-1604D - AM 
-4355LM-MG 



-^02 - MS 



-4626C - MA 
-4605 - AM 



-1355SLO-MG 
_1091H-MG 

—1438 - SE 
— 1089E - MG 
— 1438H-SE 
— 701B-SP 
-1648 - SP 
-1646 - SP 
-566 - PR 
-1357E - MG 
-1641- PR 



-fOlA-SP 
_HD29 Btgalleriare 



-USA 



-T23 Bt japonensis- Franca 
-HDl Bt kurstaki - USA 
-HD73 Bt kurstaki - USA 
_HD11 Btaizawai - USA 
-J07 Bt aizawai - Fran9a 
348L Bt galleriae - PR 
348B Btalesti-VK 
HD7 Bt dendrolimus -USA 
T16 6/ Indiana - Franca 
474 Bt galleriae -PR 
-462A Bt galleriae - PR 
""460 Bt darmstadiensis - PR 
-"461A Bt tolworthi - PR 
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~HD2 Bt thuringiensis - USA 
-1119C-G0 
-1139K~GO 
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-1132E-GO 
J135B - GO 
"Ll 148F - GO 
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J— 1131A-GO 
Lrll32A-GO 

■1109N-GO 
r'-1132C-G0 
Lrll45B - GO 
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Figure 1 Dendrogram of genetic similarity of Bacilllus thuringiensis serovars based on AFLP analysis using unweighted pair group 
arithmetic mean (UPGMA) program for 5 primer combination. 
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diversity and genetic relationship within the species of 
tropical isolates of B. thuringiensis. A total of 495 scor- 
able fragments were generated, ranging from 50 to 
500 bp, in the 65 B. thuringiensis strains when 5 primer 
combinations were used. Out of 495 fragments 483 
were found to be polymorphic, and only 12 fragments 
were monomorphic (Table 2), and 7 of these fragments 
were generated by the primer combination EcoRl-G + 
Msel-CC Our results showed that the £coRI primer 
with selective nucleotide Msel-C and Msel-CG gener- 
ated the greatest number of fragments, 130 and 113, 
respectively. However, £coRI primer with selective nu- 
cleotide AC and C with selective nucleotide Msel-CG 
generated fewer fragments, 77 and 79, respectively. 

Genetic distances were obtained using the comple- 
ment of the Jaccard similarity coefficient (CSJ), and the 
grouping using hierarchical UPGMA method. Genomic 
AFLP markers were analyzed through UPGMA cluster 
analysis to determine the genetic relationship among 
B. thuringiensis isolates (Figure 1). Out of 26 subspecies, 
strains 344 and T09 (B. thuringiensis subsp. tolworthi) 
showed the highest similarity (15%), while isolates HD3 
B. thuringiensis subsp finitimus and T24 B. thuringiensis 
subsp neoleonensis were the most genetically distant 
(92%). B. thuringiensis isolates with no subspecies identifi- 
cation, found in samples from Goias State showed higher 
similarity forming a group with an average distance of 6%, 
and the closest subspecies to this group was B. thuringien- 
sis subsp thuringiensis (HD2) with 52% of similarity. 
Groups with high similarity were formed with B. thruin- 
giensis strains collected in the same State, and other 
groups composed with B. thruingiensis strains collected in 
different Brazilian regions. No polymorphism was found 
as a reliable marker among Bt strains efficient and not effi- 
cient against fall armyworm. 

Among the 65 B. thuringiensis strains analyzed, 15 
formed a single major cluster, and coincidentiy these strains 
were sampled in the same Brazilian State, however more 
than 200 km apart. The phylogenetic analysis of these 
BrazOian soil isolates (Figure 1) illustrates the great genetic 
diversity among strains in the same serovar. However, 
B. anthracis in a previous analysis of a global collection of 
78 strains showed little variation among the isolates. 

Kumar et al. 1996 found that the low level of similarity 
between different B. thuringiensis subspecies or strains 
used in his research implied in high levels of genetic 
diversity which was reflected by the diversity in flagellar 
H-antigen agglutination reactions and presence of differ- 
ent toxins with different insect specificities. According 
to Aronson et al. 1986, several factors may contribute to 
the genetic diversity of B. thuringiensis strains, such as 
the presence of many different plasmids in size and 
numbers in each strain, and transfer of plasmids via cell 
mating. Fagundes et al. 2011, found that the number and 



size of plasmids varied between different strains in the 
same serovar, however the isolates from Goi'as State 
showed similar plasmid profiles, and these informations 
suggest that these isolates harbor common genes. Carlson, 
et al. 1996 also reported that the presence or absence of 
one or more large plasmids and they may be linear or 
circular during analysis. B. anthracis harbors two plas- 
mids and plasmids of similar size have been reported in 
B. cereus and B. thuringiensis. These observations sup- 
port the idea that horizontal gene transfer of plasmid 
and/or other extrachromosomal markers is an import- 
ant factor in defining the phenotypes of type I bacillus 
B. cereus-like isolates evolved along apparent large evo- 
lutionary distances to give rise to clusters that in more 
recent times acquired plasmids that conferred insecti- 
cidal or other pathogenic phenotypes (Hill et al. 2004). 
These observations of genetic diversity were also 
reflected in our results. Lereclus et al. 1984 mentioned 
that the transposon-like inverted repeats flanking the 
endotoxin genes are also a factor of genetic diversity. 

Few studies have been published using the fAFLP tech- 
nique to characterize B. thuringiensis isolates, but we con- 
cluded that AFLP is a highly reproducible, easy to use and 
relatively fast method which can be applied to different 
organisms, which can be used for genetic diversity studies, 
and seems particularly well suited for screening a large 
numbers of isolates. Moreover, the same protocol can be 
used in screening at the subspecies level, and demon- 
strated the value of AFLP technique for the study of diver- 
sity between strains of B. thuringiensis. 
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